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Abstract 

The  lithium /iodine  battery  is  the  most  widely  used  power  source  for  implantable  cardiac  pacemakers.  While  the  average  power  demands 
of  pacemakers  are  very  low,  the  instantaneous  power  requirements  can  challenge  the  capabilities  of  this  high  impedance  battery.  Thus,  there 
is  a  need  to  predict  the  transient-response  behavior  of  new  battery  designs.  Impedance  spectroscopy  was  used  to  characterize  the  behavior  of 
a  variety  of  cell  designs  as  a  function  of  cell  design,  depth-of-discharge,  and  current  density.  An  equivalent-circuit  model  was  developed  in 
terms  of  the  physical  elements  which  contribute  to  electrical  resistance,  capacitance,  and  electrochemical  polarization.  The  transient-response 
predicted  by  this  simple  equivalent-circuit  model  was  compared  to  pulse  data  from  actual  batteries  and  found  to  be  in  excellent  agreement. 
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1.  Introduction 

The  lithium/ iodine  battery  has  been  used  as  a  power  source 
for  implantable  cardiac  pacemakers  for  over  twenty  years.  Its 
high  reliability  and  low  self-discharge  rate  make  it  ideally 
suited  for  this  application.  As  a  result,  virtually  all  cardiac 
pacemakers  today  are  powered  by  the  lithium/ iodine  battery. 

The  lithium/iodine  battery  is  also  characterized  by  very 
high  internal  impedance  which  increases  throughout  dis¬ 
charge.  Its  characteristic  knee-shaped  discharge  curve  results 
primarily  from  the  corresponding  IR  losses.  In  the  past,  the 
power  requirements  for  cardiac  pacemakers  were  very  low 
(around  50  /zW)  and  the  high  internal  impedance  of  the  cells 
was  not  a  significant  factor  in  determining  the  useful  life  of 
the  battery.  Nearly  all  of  the  available  iodine  was  consumed 
by  the  time  the  load  voltage  dropped  to  unacceptable  levels. 

While  the  average  power  requirements  of  cardiac  pace¬ 
makers  have  tended  to  decrease,  the  instantaneous  power 
demands  have  grown  substantially  due  to  the  incorporation 
of  more  complex  electronic  circuits,  new  features,  and  new 
therapies.  As  the  instantaneous  power  demands  increase,  it  is 
possible  to  drive  the  battery  voltage  prematurely  to  unac¬ 
ceptably  low  levels.  If  the  transient-response  characteristics 
of  the  battery  are  well  understood,  many  of  these  high-power 
features  can  be  implemented  through  careful  power  manage- 
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ment.  In  general,  this  requires  that  the  battery  be  described  in 
terms  of  an  equivalent-circuit  model  and  that  such  a  circuit 
can  be  predicted  as  a  function  of  discharge  depth  prior  to  the 
development  of  a  new  battery  design. 

Previously  we  have  described  in  general  terms  the  devel¬ 
opment  of  three  types  of  models  used  to  characterize  implant¬ 
able  medical  batteries  [  1  ] .  We  have  described  in  detail  the 
development  of  a  model  which  allows  the  discharge  charac¬ 
teristics  of  a  new  cell  design  to  be  accurately  predicted  [2- 
4].  This  model  relates  the  physical  elements  of  the  battery 
design  to  the  electrical  resistance  and  electrochemical  polar¬ 
ization  observed  during  constant-current  discharge.  It  was  the 
goal  of  the  present  work  to  extend  the  model  to  predict  tran¬ 
sient-response  behavior. 


2.  Background 

For  the  sake  of  continuity  and  clarity  a  brief  review  of  the 
major  sources  of  voltage  loss  in  the  lithium  /iodine  battery  is 
repeated  here.  Details  of  the  lithium/iodine  battery  design 
and  other  characteristics  can  be  found  elsewhere  [5-7] . 

The  voltage  losses  observed  during  discharge  of  a  lithium/ 
iodine  battery  can  be  conveniently  categorized  as  ohmic 
losses,  non-ohmic  polarization,  and  loss  of  open-circuit 
potential.  The  ohmic  (current  independent)  losses  are  dom¬ 
inated  by  the  electrolyte  resistance  (i.e.,  lithium  iodide  for- 
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Fig.  1 .  Typical  impedance  spectrum  of  a  lithium/ iodine  battery  presented  in  the  complex-plane. 


mation)  during  the  first  70%  to  80%  of  discharge.  The 
cathode,  which  is  a  two-phase  material  consisting  of  crystal¬ 
line  iodine  and  an  iodine /poly-2- vinyl  pyridine  complex, 
provides  a  relatively  minor  contribution  to  the  ohmic  resis¬ 
tance  as  long  as  the  crystalline  iodine  is  available  to  saturate 
the  liquid  phase.  After  the  crystalline  iodine  is  depleted,  the 
cathode  resistance  grows  exponentially  and,  depending  on 
cell  design,  may  eventually  dominate  the  total  ohmic 
resistance. 

The  non-ohmic  (current  dependent)  resistance  arises  pri¬ 
marily  at  the  cathode /current  collector  interface  [2].  It 
should  be  noted  that  the  stainless  steel  battery  case  serves  as 
the  current  collector  in  our  cell  designs.  The  non-ohmic  resis¬ 
tance  measured  by  impedance  spectroscopy  exhibits  a  current 
dependence  consistent  with  a  charge  transfer  resistance.  In 
addition,  a  part  of  the  non-ohmic  polarization  is  related  to 
concentration  polarization  within  the  cathode.  When  the  cath¬ 
ode  composition  is  in  the  two-phase  region  (i.e.,  crystalline 
iodine  is  still  available),  the  non-ohmic  polarization  remains 
relatively  constant.  As  the  cathode  becomes  depleted  of 
iodine,  the  non-ohmic  polarization  increases  rapidly  and  is 
comparable  in  magnitude  to  the  ohmic  losses. 

The  open-circuit  voltage  of  the  lithium/iodine  battery  is 
constant  as  long  as  the  cathode  composition  is  in  the  two- 
phase  region.  After  depletion  of  the  crystalline  iodine,  the 
open-circuit  voltage  decreases  as  the  cathode  is  further 
discharged.  The  total  loss  of  open-circuit  voltage  is  typically 
less  than  30  mV  at  the  point  at  which  the  load  voltage  falls 
to  unacceptable  levels.  Thus  the  loss  of  open-circuit  voltage 
represents  the  least  significant  source  of  voltage  loss  in  a 
practical  lithium /iodine  battery. 

Fig.  1  shows  a  typical  impedance  spectrum  of  a  lithium/ 
iodine  battery  in  which  the  ohmic  and  non-ohmic  contribu¬ 
tions  to  the  cell  impedance  are  noted. 


3.  Experimental 

The  batteries  used  in  these  studies  were  hermetically 
sealed,  central  anode  designs  [5].  Nine  separate  battery 
designs  were  studied  having  capacities  ranging  from  about 
0.60  to  about  1 .7  Ah.  Discharge  current  densities  ranged  from 
about  2  to  about  10  mA  cm  2.  Impedance  spectra  were  col¬ 
lected  with  a  Solartron  1255  frequency  response  analyzer 
interfaced  to  an  EG&G  PAR  273a  potentiostat.  Typically, 
the  cells  remained  on  load  during  the  measurement  of  impe¬ 
dance  data. 


4.  Results  and  discussion 

As  is  typical  of  electrochemical  systems,  the  complex- 
plane  plots  contain  somewhat  flattened  semi-circles  (see Fig. 
1 )  and  cannot  be  perfectly  represented  by  a  pure  resistor  and 
capacitor  in  parallel.  For  our  purposes,  however,  it  was  nec¬ 
essary  to  approximate  the  impedance  behavior  using  ideal 
circuit  elements.  This  allows  the  battery  behaviour  to  be  sim¬ 
ulated  through  commercially  available  circuit  simulation 
software  or  to  be  approximated  with  the  use  of  a  power  supply 
and  a  collection  of  resistors  and  capacitors.  These  approaches 
allow  system  simulations  and  bench  tests  to  be  completed 
prior  to  finalization  of  circuit  and  battery  design.  Our 
approach  was  to  use  the  intersections  of  an  arc  with  the  real 
axis  to  define  the  resistance  value  and  to  calculate  the  capac¬ 
itance  as: 

1 


where  wmax  is  the  angular  frequency  at  the  apex  of  the  arc. 
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Fig.  2.  Double-layer  capacitance  (normalized  for  anode  area)  as  a  function  of  discharged  capacity  for  a  single  battery  design. 


The  effect  of  approximating  the  impedance  behavior  in 
this  manner  will  be  discussed  later. 

4. 1.  Double-layer  capacitance 

As  mentioned  in  Section  2,  the  non-ohmic  impedance  has 
been  shown  to  be  consistent  with  a  charge- transfer  resistance 
residing  at  the  cathode/case  interface.  Thus  it  would  be 
expected  that  the  corresponding  capacitance  represents  a  dou¬ 
ble-layer  capacitance.  Furthermore,  it  might  be  expected  that 
the  double-layer  capacitance  would  be  relatively  independent 
of  cathode  composition  and,  hence,  independent  of  discharge 
depth. 

Fig.  2  shows  the  measured  capacitance  associated  with  the 
charge  transfer  resistance  for  a  given  cell  design  as  a  function 
of  depth  of  discharge.  These  data  represent  approximately  40 
cells  discharged  at  current  densities  ranging  from  2  to  5  mA 
cm  “  2.  The  measured  capacitance  is  normalized  by  the  anode 
area  under  the  assumption  that  the  close  proximity  of  the 
anode  to  the  case  would  tend  to  localize  the  current  density 
to  this  area. 

The  average  value  of  the  area  normalized  capacitance  is 
around  10  ^uF  cm-2.  This  is  a  reasonable  value  for  a  double¬ 
layer  capacitance  and  provides  support  for  this  hypothesis. 
Although  it  is  not  apparent  from  the  data  as  presented  in  Fig. 
2,  there  is  no  correlation  between  the  measured  capacitance 
and  the  discharge  rate.  The  scatter  in  the  data  represents 
normal  cell-to-cell  variation.  There  does,  however,  appear  to 
be  some  decrease  in  the  measured  capacitance  at  greater 
depths  of  discharge.  A  similar  trend  was  observed  for  every 
cell  design  included  in  this  study,  except  that  the  cells  with 
less  capacity  exhibited  steeper  slopes. 

Fig.  3(a)-(d)  shows  the  double-layer  capacitance  of  sev¬ 
eral  cell  designs  plotted  as  a  function  of  depth-of-discharge 


instead  of  capacity.  These  cells  were  of  a  similar  design  in 
that  they  were  the  same  thickness  and  width  and  differed  only 
in  the  height  dimension.  The  anode  area  and  capacity  of  each 
cell  is  given  in  the  legend.  The  capacitance  was  again  nor¬ 
malized  by  the  anode  area.  As  in  Fig.  2,  there  is  significant 
scatter  in  the  data  representing  typical  cell-to-cell  variation. 
However,  each  plot  yields  similar  area  normalized  capaci¬ 
tance  values  and  a  similar  downward  trend  in  capacitance 
with  depth-of-discharge.  Nearly  identical  results  were 
obtained  with  all  cell  designs  included  in  this  study. 

One  possible  explanation  for  this  behavior  may  be  that  the 
normalization  of  the  double-layer  capacitance  by  the  anode 
area  may  not  be  strictly  correct.  Because  the  current  densities 
are  low  and  the  interfacial  resistance  is  high,  the  effective 
area  of  the  cathode/  case  interface  may  be  significantly  greater 
than  the  area  of  the  adjacent  anode.  The  total  cathode/case 
interfacial  area  exceeds  the  anode  area  by  a  factor  of  1 .7  to 
2.4  depending  on  cell  design.  As  the  cathode  is  discharged, 
its  volume  decreases  and  the  remaining  cathode  becomes 
located  immediately  over  the  anode.  By  the  time  discharge 
of  the  cathode  is  nearly  complete,  the  area  around  the  edges 
of  the  anode  are  completely  devoid  of  cathode.  Thus  it  may 
be  more  appropriate  to  normalize  the  double-layer  capaci¬ 
tance  by  the  entire  cathode/ case  interfacial  area  at  the  begin¬ 
ning  of  discharge  and  by  the  anode  area  at  the  end  of 
discharge.  This  would  result  in  an  decrease  in  capacitance  by 
a  factor  of  1.7  to  2.4  between  the  beginning  and  end  of 
discharge.  This  agrees  approximately  with  the  observation  in 
Fig.  3(a)-(d). 

4.2.  Geometric  capacitance 

The  remaining  capacitance  can  be  assigned  to  the  geomet¬ 
ric  capacitance  between  the  central  lithium  anode  and  the 
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Dopth  of  Discharge  Depth  of  Discharge 

Fig.  3.  Double-layer  capacitance  (normalized  for  anode  area)  as  a  function  of  depth-of-discharge  for  several  battery  designs  differing  in  one  dimension.  Anode 
areas  and  cathode  capacities  are:  (a)  3.8  cm2, 610  mAh;  (b)  5.6  cm2,  830  mAh;  (c)  7.3  cm2,  1000  mAh;  (d)  10.8  cm2,  1300  mAh. 

battery  case.  Fig.  4  shows  how  the  geometric  capacitance  There  are  at  least  two  factors  which  contribute  to  the 

varies  throughout  discharge  for  a  typical  cell  design.  Note  decreasing  value  of  the  geometric  capacitance.  First,  as  the 

that  the  magnitude  of  the  geometric  capacitance  decreases  by  lithium  anode  is  consumed,  the  anode-to-case  spacing 

about  two  orders  of  magnitude  from  the  beginning  to  end  of  increases  by  a  factor  of  about  three.  Second,  the  composition 

discharge.  of  the  dielectric  material  changes  radically.  At  the  beginning 
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Fig.  4.  Geometric  capacitance  as  a  function  of  discharged  capacity  for  a  typical  cell  design. 
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Fig.  5.  Charge-transfer  resistance  and  total  non-ohmic  polarization  as  a  function  of  discharged  capacity  for  a  typical  cell  design. 


of  discharge,  the  dielectric  material  consists  of  a  very  thin 
film  of  lithium  iodide  and  a  much  thicker  layer  of  the  cathode 
material. 

By  the  end  of  discharge,  the  inter-electrode  space  is  filled 
almost  entirely  by  lithium  iodide.  While  we  have  not  meas¬ 
ured  the  dielectric  constant  of  the  cathode  material,  we  expect 
it  to  be  much  higher  than  that  of  lithium  iodide. 

We  have  collected  much  less  data  for  the  geometric  capac¬ 
itance  than  for  the  double-layer  capacitance.  One  reason  for 
this  is  that  the  time  constant  for  the  geometric  capacitance 
and  ohmic  resistance  is  small  enough  that  the  effect  of  the 
geometric  capacitance  only  becomes  relevant  at  frequencies 
above  about  10  000  Hz.  In  most  of  our  applications,  we  are 
not  interested  in  the  transient-response  of  the  battery  at  these 
high  frequencies.  The  second  reason  for  our  limited  data  in 
this  area  is  that  the  collection  of  accurate  high  frequency  data 
requires  extreme  care  and  precludes  the  use  of  a  potentiostat. 
This  process  is  cumbersome  and  has  been  avoided  wherever 
possible. 

Our  approach  for  modelling  the  geometric  capacitance, 
when  needed,  has  been  to  fit  an  empirical  function  to  the 
limited  data  and  then  scale  the  data  by  the  anode  area. 

4.3.  Ohmic  resistance 

As  described  earlier,  the  ohmic  resistance  is  determined  by 
the  electrolyte  resistance  and  cathode  resistance.  It  is  not 
possible  to  resolve  these  two  components  by  impedance  spec¬ 
troscopy  at  any  depth-of-discharge.  A  description  of  the  func¬ 
tional  forms  used  to  model  the  ohmic  resistance  can  be  found 
in  Ref.  [4]  and  will  not  be  discussed  further  here. 


4.4.  Non-ohmic  polarization 


In  our  previous  work,  we  modelled  the  entire  non-ohmic 
polarization  with  a  single  function  which  incorporated  loss 
of  open-circuit  potential,  concentration  polarization,  and  acti¬ 
vation  overpotential.  This  was  a  useful  approach  for  predict¬ 
ing  the  load  voltage  of  the  battery  under  constant  current 
discharge  conditions.  It  is  not  useful,  however,  for  predicting 
transient-response  characteristics  since  each  of  these  contri¬ 
butions  has  a  drastically  different  time  dependence. 

As  discussed  earlier,  the  open-circuit  voltage  is  constant 
throughout  most  of  discharge  and  amounts  to  no  more  than 
about  30  mV  at  the  end  of  discharge.  The  remainder  of  the 
non-ohmic  polarization  is  defined  by: 


^7conc 


T  ^7act  Voc  3R 


ohmic 


or  by: 


L  non-ohmic 


ohmic 


Both  the  ohmic  resistance,  Rohmic,  and  the  charge- transfer 
resistance,  Rct ,  are  measurable  by  impedance  spectroscopy  as 
shown  in  Fig.  1 .  The  d.c.  resistance,  Rdc,  is  calculated  from 
the  load  voltage  and  the  concentration  polarization,  expressed 
as  either  7]conc  or  Rconc,  is  calculated  by  difference. 

Fig.  5  compares  the  charge  transfer  resistance  to  the  total 
non-polarization  as  a  function  of  discharged  capacity  for  a 
typical  cell  design.  As  expected  from  our  previous  work,  both 
quantities  remain  relatively  constant  until  the  cathode  reaches 
its  single-phase  composition.  Then  both  quantities  increase 
rapidly  as  the  iodine  activity  decreases. 

The  most  surprising  result  of  this  study  is  shown  in  Fig.  6. 
Here  the  ratio  Rct/ ^non-ohmic  is  plotted  versus  discharged 
capacity  for  a  typical  cell  design.  This  particular  plot  repre- 
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Fig.  6.  Ratio  of  charge-transfer  resistance  to  total  non-ohmic  polarization  as  a  function  of  discharged  capacity  for  a  typical  cell  design. 


sents  the  same  set  of  batteries  used  to  create  Fig.  2.  Note  that 
the  impedance  data  were  collected  over  a  range  of  current 
densities.  Although  there  is  a  good  deal  of  cell-to-cell  varia¬ 
tion,  the  ratio  ^Cc  ^non-ohmic  is  essentially  independent  of  cur¬ 
rent  and  cathode  composition.  The  average  value  of  the  ratio 
is  about  0.45.  We  have  observed  nearly  identical  behavior 
and  a  similar  ratio  for  all  of  the  cells  designs  included  in  this 
study. 

It  is  not  surprising  that  both  Rcl  and  Rconc  increase  with 
decreasing  iodine  activity  or  that  both  decrease  with  increas¬ 
ing  current  density  ( as  long  as  the  current  density  is  far  below 
its  mass  transport  limit) .  Conversely,  it  is  very  surprising  that 
the  ratio  remains  so  constant  over  such  a  wide  range  of  current 
density,  cathode  composition  and  cell  design.  Although  we 
do  not  fully  understand  this  phenomenon,  this  simple  rela¬ 
tionship  allows  a  convenient  means  to  estimate  the  charge 
transfer  resistance  from  the  total  non-ohmic  resistance  as 
calculated  from  the  function  in  Ref.  [4] . 

4.5,  Equivalent- circuit  model 

Fig.  7  shows  an  approximate  equivalent-circuit  model  for 
the  lithium/iodine  battery  which  summarizes  all  of  the  com- 


(Rel  +  Rcath) 

Fig.  7.  Equivalent  circuit  model  for  the  lithium /iodine  battery. 


ponents  described  above.  Note  that  we  have  chosen  to  rep¬ 
resent  the  concentration  polarization,  rjconc,  as  voltage  source 
opposed  to  the  open-circuit  voltage,  Voc.  The  time  constant 
for  the  development  of  concentration  polarization  is  of  the 
order  of  minutes  to  hours.  Thus,  while  it  is  a  significant 
contributor  to  the  voltage  loss  observed  during  a  constant 
current  discharge,  it  can  be  treated  as  a  time- independent 
voltage  loss  on  a  short  time  scale. 

4.6.  Comparison  of  model  to  pulse  data 

When  developing  equivalent-circuit  models  for  batteries, 
it  is  imperative  that  the  models  be  verified  for  the  particular 
range  of  applications.  In  this  case,  for  example,  it  is  necessary 
to  determine  if  the  use  of  ideal  circuit  elements  to  represent 
the  observed  impedance  behavior  is  an  acceptable  approxi¬ 
mation.  Additionally,  it  must  be  shown  that  ignoring  the  time 
dependence  of  concentration  polarization  does  not  introduce 
a  significant  error  on  the  time  scale  of  interest.  One  compar¬ 
ison  of  the  model  to  impedance  and  pulse  data  is  shown  in 
Figs.  8-10. 

Fig.  8  compares  a  complex-plane  plot  of  a  measured  impe¬ 
dance  spectrum  to  a  corresponding  impedance  spectrum  sim¬ 
ulated  with  ideal  circuit  elements.  Fig.  9(a)  and  (b)  shows 
the  same  information  on  a  Bode  diagram.  While  the  differ¬ 
ences  in  the  complex-plane  plot  are  striking,  inspection  of  the 
Bode  diagram  shows  that  much  of  the  difference  lies  in  the 
measured  versus  calculated  phase  shift  at  certain  frequencies. 
The  corresponding  differences  in  the  magnitude  of  the  meas¬ 
ured  and  simulated  impedance  are  relatively  small. 

The  conformance  of  predicted  pulse  behavior  to  observed 
behavior  is  explored  in  Fig.  10(a),  (b)  and  (c).  These  data 
were  collected  during  the  application  of  a  50  mA  pulse  to  the 
battery.  Fig.  10(a)  shows  the  excellent  agreement  between 
the  model  and  the  observed  behavior  at  short  time  frames.  At 
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Fig.  8.  Comparison  of  measured  and  simulated  impedance  spectra  in  the 
complex-plane. 

slightly  longer  times  (Fig.  10(b)),  the  model  predicts  a 
slightly  lower  voltage  than  is  actually  observed.  This  slight 
difference  arises  from  the  use  of  an  ideal  capacitor  in  the 
equivalent-circuit  model.  At  still  longer  times  (Fig.  10(c)), 
the  observed  voltage  drops  below  the  predicted  voltage. 

This  demonstrates  the  onset  of  concentration  polarization 
which  is  not  treated  as  a  time-dependent  quantity  in  the  equiv¬ 
alent-circuit  model. 

Overall,  there  is  reasonably  good  agreement  between  the 
observed  and  predicted  behavior  over  a  time  period  of  about 


Timc/s 


Fig.  9.  Comparison  of  measured  and  simulated  impedance  magnitude  (a) 
and  phase  angle  (b). 


100  s.  At  significantly  longer  times,  the  constant-current 
model  described  in  Ref.  [4]  becomes  applicable.  Thus,  by  a 
combination  of  these  two  models  it  is  possible  to  predict  the 
load  voltage  of  new  battery  designs  as  a  function  of  discharge 
capacity  over  a  wide  range  of  constant-current  and  transient 
conditions.  Furthermore,  it  is  possible  to  estimate  perform¬ 
ance  distributions  in  the  transient  model  in  a  manner  similar 
to  that  we  have  described  for  the  constant-current  case.  The 
performance  variation  measured  in  the  cells  in  this  study 
provides  an  excellent  description  of  the  distribution  functions 


Tim*  /  * 

Fig.  10.  Comparison  of  measured  and  predicted  voltage  response  to  a  50  mA  pulse  over  the  intervals:  (a)  0-0.1  s;  (b)  0-1  s;  (c)  0-100  s. 
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for  each  of  the  equivalent-circuit  elements.  This  information 
can  be  used  in  Monte  Carlo  simulations  to  extend  the  model 
to  estimate  performance  distributions  for  new  cell  designs. 
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